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We point out the existence of a new type of growing transverse mode in the gravitational 
instability. This appears as a post-Newtonian effect to Newtonian dynamics. We demon¬ 
strate this existence by formulating the Lagrangian perturbation theory in the framework of 
the cosmological post-Newtonian approximation in general relativity. Such post-Newtonian 
order effects might produce characteristic appearances of large-scale structure formation, 
for example, through the observation of anisotropy of the cosmic microwave background 
radiation (CMB). 


¥*) 


§1. Introduction 


The investigation of nonlinear, large-scale structure formation of the universe 
has been a major research subject in cosmology. As long as the scales of the density 
fluctuations are much smaller than the horizon scale, the Newtonian approximation 
is sufficiently accurate to describe their evolution up to the non-linear regime, where 
the density contrast becomes larger than unity. On the other hand, the region in 
which the large-scale structure is observed is steadily growing. For example, SDSS 
(Sloan Digital Sky Survey) will cover over a region of several hundred megaparsecs. 
It is not clear if the application of Newtonian theory is appropriate for such a wide 
region of spacetime. In fact fluctuations relevant for such large-scale structures are 
not always much smaller than horizon scales in the past. For example, the horizon 
scale at the decoupling time is on the order of ctdec(l + -z^dec) ~ 80/i“^Mpc in the 
present physical length. This suggests that we must employ a relativistic description 
for the evolution of fluctuations larger than or equivalent to such scales. 

Thus to understand the evolution of the large-scale structure of the universe, 
it is important and necessary to have some formalism to clearly evaluate the effect 
of general relativistic corrections to the Newtonian dynamics. Naturally, we also 
require the formalism to agree with the gauge invariant linear theory developed by 
Bardeen&,and Kodama and SasakiS^ in the linear regime. 

We develop such a formalism based on the Lagrangian perturbation theory in 
the post-Newtonian (PN) framework in which it has been shown that the formal¬ 
ism can be app^d even for perturbations larger than the present horizon scale 
~ 2000/i“^Mnc.EiI The Newtonian Lagrangian picture has also been developed by 
Buchert,l3)’Q)’l3 where it is expected that the aapmaiih gives a good approximation 
up to a certain stage of the non-linear regime. There have also been some 
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studies based on relativistic Lagrangian perturbation theory. However, because 

of the gauge condition adopted in these studies, namely, the synchronous comoving 
coordinates, it is not easy to have contact with the Newtonian Lagrangian approach 
fully developed and used for the numerical simulation. Thus we shall study the La¬ 
grangian perturbation in coordinates where comparison with the Newtonian case is 
most easily carried out. 

§2. Cosmological post-Newtonian approximation 


As stated above, we wish to develop the Lagrangian perturbation theory formu¬ 
lated in a coordinate system which has a straightforward Newtonian limit. For this 
purpose it is convenient to use the (3-1-1) formalism motivated by the following con¬ 
siderations. Let us first assume that there exists a congruence of time-like worldlines 
from which the spacetime looks isotropic. We shall call the family of the worldlines 
basic observers, who see no dipole component of the cosmic microwave background 
radiation (CMB). We can regard that any one of these observers at the spacetime 
point X moves with 4-velocity n^(x), without loss of generality. The tangent vector 
is normalized as = — 1. These observers are used to foliate the spacetime 

by their simultaneous surfaces: t = const. These considerations imply that matter, 
such as a galaxy, moves with a velocity in our coordinates according to the geodesic 
equation. Hence we naturally expect that the dynamics of nonrelativistic matter is 
described by the cosmological Newtonian picture at the lowest order. Fortunately, 
there have been some studies on the post-Newtonian approximation applied to cos¬ 
mology©’ 0 in the (3-1-1) formalism. We shall follow these approaches to work in 
the following coordinates: 

ds^ = —(a^ — PiP^){cdt)‘^ -|- 2Pi{cdt)dx^ -|- a^(t)[(l — -|- hijjdx'^dx^, (1) 


where = 0. We employ gauge conditions termed the “longitudinal gaugi 

or “Newtonian gauge” with regard to the scalar modes of the metric perturbations: 


= 0 , = 0 . 


( 2 ) 


Here ' denotes the covariant derivative with respect to the background met- 
ric yb , for which we take the spatial metric of the Friedmann-Robertson-Walker 
(FRW) geometry as = 6ij/{l + /Cr^/4)^. The quantity /C is the curvature pa¬ 
rameter of FRW models and = xf + X 2 + x'^. We treat the metric perturbations 
as small quantities, so we are able to regard tensorial quantities in our equations 
as tensors with respect to the background spatial metric yb . The above condition 
(2) guarantees that hij and Pi contain only the tensor mode, which represents the 
freedom of the gravitational wave in the PN order, and the vector mode, respectively. 

The explicit forms of^iuatein equations in terms of the above variables can be 
found in previous works. the following, we only consider the Einstein-de 

Sitter background universe (/C = A = 0) and use Cartesian coordinates for the sake 
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of simplicity: 

Hit? = 

\a{t)) 3 


(3) 


where H{t) and pb{t) are the Hubble parameter and the homogeneous density of the 
background FRW universe, respectively. For simplicity, we assume that the above 
equation is derived by the usual averaging method 0^ and we may regard pb as an 
average value over the volume as large as the horizon scale (ct)^: {p)(^ctp = Ph{t)- 
The metric required for the usual Eulerian Newtonian picture in the perturbed 
FRW universe is known to take the following form: 


ds"^ 



c^d? + a^{t) ^1 



Sijdx^'dx^. 


(4) 


The quantity (pN is a cosmological Newton-like gravitational potential related to the 
matter density fluctuation field <5 at via the Poisson equation 


Ax4>n = 


(5) 


where SNix, t) = {pn{x, t) — Pb{t))/Pb{t) , and is the Laplacian. The above metric 
is usually used to give an accurate description of the trajectories of nonrelativistic 
fluid elements on scales much smaller than the Hubble distance cH~^. 

We expand the basic equations in powers of c~^ to obtain the post-Newtonian 
approximation with the condition that the lowest order metric takes the above New¬ 
tonian form (Q). Thus the PN terms of all metric variables used in this paper may 
be expanded as follows: Edl) 


a = 1 + 4 + ^ + 0(c-6), 




■ d®* 

/9* = ^ + 0(c-'), 


h. 


(4) 

d 


hij = -h 0(c ^), 


( 6 ) 


where cj) is the Newtonian gravitational potential as stated below. It should be noted 
that we consider only the PN expansion of /3*, not of /3j. If we adopt a dust model 
= p?Uf^u,^ for the energy-momentum tensor of matter, the metric up to order c~^ 
agrees with (Q) by substituting the above expressions into Einstein equations. Erom 
the other substitutions, we can derive the relations between the metric perturbations 
and the matter variables in each order of c“”. We shall here present only relevant 
equations in our calculations. The lowest order field equation gives us 


A^fj) = A^'ip^‘^'> = A-irGa^ip - pb). 


(7) 


The field equations at the next order, that is, at the first PN order take the following 
forms for and 


Axa 


(4) 


4),k (p,k +47rGa^ (2pa‘^v‘^ 


pcf) -I- 3/96(/>^ - 3a^ 


d? a dt J 


( 8 ) 
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= leyrGa^pu* + 4 + ^4>,i ^ , 


(9) 


where we have used the Newtonian order equations (7) and the gauge conditions (2), 
and n® is defined as u®/c = u^ju^. We note that u* represents a peculiar velocity of 
a fluid element in the comoving frame, so the physical peculiar velocity is an®. 
Similarly, the material equations up to the first PN order become 


dt 


pa-^ {1 + ^ ~ 


+ 


A 

9x® 


pa'^v"- ^ 1 + ^ ~ 


2 2 
' V 

^-4 


+O(c-^)=0, (10) 


(9n® ■ 9n® d ^ 1 ^ 

+ + 2-n® + ^n® 






a /I 


dt \2 


— 3(/> + n-^ 


a /I 


— 3(^ 


-A .-11 

9 9 9 


dx^ \2 
d 


a^^\i +3())(?i,i +a^v^(l),i +^(a^ A®) 

dt 




+ 0(c-^).(ll) 


The lowest order in these two equations reduces to the Newtonian Eulerian equations 
of hydrodynamics for a pressureless fluid, and the terms of order c~^ provide the first 
PN corrections. Equations (10) and (11) are our basic equations for the Lagrangian 
approach to the trajectory field of matter fluid elements. 


§3. PN Lagrangian perturbation approach to large-scale structure 
formation in the perturbed FRW universe 


3.1. Basic equations for the transverse part of the trajectory field of the fluid ele¬ 
ments 

We shall rewrite the above set of post-Newtonian equations in the Lagrangian 
coordinates and then solve them perturbalively by extending the formalism in the 
Newtonian theory developed by Buchert.il® Namely we concentrate on the integral 
curves x = f{X,t) of the velocity field v{x,t): 


A 

dt 




v{f,t), 


f{X,ti) = X, (12) 


where X denotes the Lagrangian coordinates which label each fluid element, x is the 
position of such an element in Eulerian space at time t, and tj is the initial time when 
the Lagrangian coordinates are defined. It should be noted that we introduce the 
Lagrangian coordinates on the comoving coordinate because we have already derived 
basic equations by using the ERW metric defined with the comoving coordinates. 
As long as the mapping f is invertible, we can give the inverse of the deformation 
tensor fl^j, which is written in terms of variables {X,t): 


dA® 

dx^ 


= hij{X,t) 


2j^ia,b^jcdfc\afd\b: 


t 


(13) 
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where J is the determinant of the deformation tensor /j|j, and is a totally 
antisymmetric quantity with £123 = +1- The comma and the vertical slash in the 
subscript denote partial differentiation with respect to the Eulerian coordinates and 
the Lagrangian coordinates, respectively. 

The continuity equation (10) then may be rewritten as 


A 

dt 


pa^ {l + ^A{X,t)) JiX,t) 


+ 0(c-4) = 0, 


(14) 


where 

AiX,t)^^a\t)f-3<PiX,t). (15) 

Thus the density field is integrated exactly up to the first PN order in the Lagrangian 
picture as in the Newtonian case: 


p{X, t)a^J{X, t) = (^1 + ^A{X,t)^ ^ P{X)a^ (^1 + ^ A (X)) + ©(c"^), 

(16) 

where the quantities with such as a, are the quantities at the initial time tj, 
Ax (f>= 47rG d {p — p^) with Ax representing the Laplacian with respect to the 

o 

Lagrangian coordinate, and J= 1. 

The trajectory field f is determined by solving the PN equations of motion 
(11). As in the Newtonian case,ll2P we consider the rotation and the divergence of 
the equations of motion with respect to the Eulerian coordinates, respectively. In the 
Lagrangian picture, the former set becomes basically the evolution equation of the 
transverse part of the trajectory field, and the latter becomes the evolution equation 
of its longitudinal part. We shall here concentrate on the transver^ part at the PN 
order. The longitudinal part will be discussed elsewhere in detail. 

So by operating with rotation on Eq. (11) with respect to the Eulerian coordinate 
we obtain 


^ijk 


d 

dx^ 


/" + 2 -/" + 4 

a 

d 


f^A 


1 1 


d 


(? a? dx^ - 


2 2 i 
a V (j)h 




According to the usual procedure in the Lagrangian formalism, we can modify the 
above equation by using Eqs. (^) and (p^ in the following form: 


^abcfj\afi\b 


fj\c + ‘^~fj\c 


J_ 

72 




^ + jeabcfi\bflfl\c4'\a 


dx^ 


2^o,bcfj\afi\b(^fj\c-^N TO^C 


(18) 


where we have not transformed the terms including the shift vector in Eq. (18) 
into their forms in the Lagrangian picture, because the expression in the Eulerian 
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picture will make it easier to understand the expansion of the above equation in later 
discussion. We note that the right-hand side of Eq. (18) is of order c“^; namely, it 
has only quantities of first PN order. 

The structure of the above equation allows us to solve for the displacement vector 
iteratively up to any desired order in in terms of the lower order displacement 
vectors and metric variables. 


3.2. The solution of the transverse part at PN order 

The solution of Eq. i) in an Einstein-de Sitter background can be chosen as 


a{t) = 



(19) 


where to denotes the present time. We use this normalization for the scale factor 
since it makes the physical interpretation of the solution more transparent, partly 
because the length in the comoving frame is equal to the physical length at the 
present universe with this choice. Eurthermore, we shall assume that the initial 
hypersurface exists at sufficiently early time so that the initial density contrast field 
is much smaller than unity. This allows us to adopt the order of the initial density 
contrast, say A, as a new perturbation parameter. 

The parameter A is assumed to be small and dimensionless. We formally split 
the initial density field accordingly: 


P(X) =Pb+APU(X), (^p{X)'j.=Pb, (20) 


where 6 denotes the initial density contrast field. This choice of the initial data is 
adequate, because the density field at some given time need not be perturbed in the 
Lagrangian framework. 

We assume that the vorticity of the Newtonian displacement vector is negligi¬ 
ble, which seems to be reasonable. We consider the following ansatz for the 

transverse perturbations at the first PN order superposed on the solved Newtonian 
irrotational trajectory field: 


f{X,t) = ^^Y{X,t) {XN)) , 

( 21 ) 

where 

^^^nx, t) = X^ + Xq,{t)F^i\^X) + X\Ut)^^Nl (22a) 


with 



qzz = 



3 

14 



/ t \ 2/3 1 4 

\Ti) 



(22b) 

(22c) 
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and 


= -I 6g, 




(1) WrP) 




N\ii^N\jj ^N\ij^N\ij- 


(22d) 

(22e) 


The initial conditions are imposed on the time coefficients by the definition of 
the Lagrangian coordinate (|l^): Q"^{tj) = — 0- quantity §g is defined as 

S —2 (/) denotes the initial gauge invariant densi^ contrast field. ©’0) N)p 

is the Newtonian solution up to the second order in A,® and we have here ignored 
the PN part of the longitudinal part (discussed elsewhere©). The PN transverse 
displacement vector at each order is constrained by 


Vx • = Vx • = 0. 


(23) 


Before considering the PN rotational flow in our coordinates, we remark that we can 
express the Newtonian gravitational potential constrained by Eq. (7) in terms of the 
above Newtonian solution following the Lagrangian perturbation formalism: 


(j)N{X,t) = -A 


tj{i + ziy ^ 


- A^ 


^2ad(jzz T O' Qzz^ 


+Qz {^odqz + 0^^(X) + O(A^), (24) 


where 


<HX) = - Uy 


47r|X-l"| 


(25) 


In the above derivation, we have used the Green function —l/(47r|Ar — 1^1) of the 
Laplacian motivated by the following considerations. As seen in Eq.(7), the Newton¬ 
like gravitational potential in the cosmological situation is generated by the density 
contrast 6. Then we can safely express such quantities using the Green function of 
the Laplacian, because the integrals over the horizon scale converge to a definite value 
by assuming that the density contrast field obeys periodic boundary condition on the 
horizon scale. We have found that the cosmological PN treatment to the large-scale 
strud^e formation is valid even for perturbations larger than the present horizon 
scale.© 

We need Eq. (18) to solve the transverse part Si of the PN displacement 
vector. Naturally, the order (P in this equation produces Newtonian counterparte in 
the case that the trajectory field is introduced on the comoving coordinates. 

Then, if we assume Jn > Jpn/c^, where Jpn is a term of order c~^ in expanding J, 
and take into account the order of A in each term on the right-hand side, Eq. (18) 
becomes 


^abcfj\afi\b 


4|c + 2-/,|c 

ji a . 

Jn 




c 

d r, 


(N)fl(fD{‘i)k _ p(3)l _ V 


I /n 


7„-4 \3j 
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where we have neglected terms like (I)nv%- in our approximation because they are 
of order at most, and Jw = det(('^Vj|j)- This equation shows that the shift 
vector constrained by the gauge condition (2) generates the first PN transverse flow. 
Inserting the ansatz (21) into Eq. (26), we obtain the following differential equation 
for solving the transverse velocity potential of the first PN displacement vector at 
the order c“^: 


A(Q^ + 2“ 


(Vx X ('’"'-I" 

Jn d 


is(i))^ + A2 fe + 2-Q,] (Vx 


{PN)„{2) 


( a.+ 2“<iA el'"'"'-''' 


l\k 




The form of this equation shows that it can be solved if we express the shift vector 
in terms of already-known Newtonian trajectory field 

Thus we first must solve the shift vector with respect to the Lagrangian coordi¬ 
nate. The equation (27) shows that we need only the explicit form of its transverse 
(3)i 

part, say (3^ , for the situation in which we are interested here. We note that 
since the gauge condition (2) implies that the transverse part of the shift vector is a 
gauge invariant quantity, no ambiguities caused by the gauge freedom remains. The 
perturbation quantity is constrained by the Einstein equation (9): 

= 4 (28) 


The right-hand side of this equation is shown to be divergenceless through the con¬ 
tinuity equation (10) at Newtonian order. We have already obtained an expression 
for the peculiar velocity field and the peculiar gravitational potential ( 0 ) and (24) 
in terms of the Lagrangian coordinate, respectively. Therefore, it is adequate to 
make use of Eq. (28) for our purpose. Then if we express Eq. (28) in terms of 
the independent variables X and t, we can obtain the following equations in simple 
form: 




= A^ 


tj{l + ziy 


-1/3 


(-^ji/t'l’i + + <11) + Ol^") . (29) 


where we have used Eqs. (|2^, (24) and the equation of mass conservation, ([l^). We 

are now in a position to be able to solve for the shift vector iteratively. Accordingly, 

(3)^ 

we can conclude that the shift vector is of second order in A and obtain the 
following equation at order A^ of Eq.(29): 


= A2 


tf(l -h Zl)"^ \tj 


-1/3 




(1) uAi) 


_i_ 1^(1) q, 

^ AJ\A ' ^ AT\AX^ AT\A ' ^ \ 


,(2) 


N\j]N\i ^ ^N\irN\j 


N\i 


(30) 
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The consistency of our formulation can easily be checked if one confirms that the 
right-hand side of this equation also satisfies the divergenceless condition with respect 
to the Lagrangian coordinate by using Eq. (22e). If we again use the Green function 
of the Laplacian, we can express the solution of this equation as 


i3P\x,t) = + o(a’) , (31) 


where 


:= d^Y 


tf(i + ziY \ti 

(r) _ 

47 r|X - Y\ 


.(32) 


In the above expression of the integration, we have again assumed that the homo¬ 
geneous solution of Eq. (30) is zero. It should be noted that we consider not 


Pi , throughout this paper. If one needs the quantity pPp one finds that it has a 
growing mode by the definition P^^'^ = The above equation is the expression 

of the shift vector that we have sought in order to solve Eq. (27). Thus we are able 
to express all perturbation quantities of the metric in our coordinates in terms of 
the trajectory field. This must be the case because it is only a dynamical variable 
in the Lagrangian description. 

By inserting Eq. (31) into Eq. (27), we can obtain the following equation at the 
lowest order of A: 


(3) 


(^Q^ + 2^Q^^Vxx(-P^)sW(X) = 0. (33) 


As an irrotational case, the form of this equation allows us to seek solutions of the 
form 


Q^ + 2-Q^ = 0, with (^^)S«(X) = iT(A:), Vx • TT = 0, (34) 

a 

where II is an unknown function determined by the initial conditions. The solutions 
of Eq. ( p^ ) can be easily found to have only a decaying mode. The first order solution 
may be safely ignored, because the decaying mode plays no physically important role: 

(JT(PA-) =.(1) ^ g (35) 

Similarly, by inserting Eqs. ( |3l| ) and (|35D into Eq. (27), we can obtain the following 
equation at second order in A: 

(qI + 2‘^Ql) Vx X <™)S« = Vx X aP. (36) 

The solution of this equation can be found by solving the linear ordinary differential 
equations 

Qr. + 2-QL = ^ , (37) 
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with Ax^^^W^ = - 


1 


(1) .7,(1) 


{l + zi)H] 


-Fir . F 


N\j3N\i 


/(I) 

' N\iF 


( 1 ) 

N\j 


+ F 


( 2 ) 

N\i 


(38) 


Since we expect the quantity to be much smaller than the order A 

o 

quantities such as So at the initial time, we can safely adopt the additional initial 
condition ^ = 0 as a longitudinal perturbation case^ besides the 

initial condition = 0. Then we can find the solution of Eq. (p^, 


Ql 




12 + 8 



(39) 


§4. Discussion and snmmary 


We have found that the transverse part of the first PN displacement vector has a 
growing mode even if the Newtonian trajectory field does not have a transverse mode. 
We emphasize that the PN solution with the growing mode is a particular solution of 
the differential equation (|^) and is generated by the initial density fluctuation field 
through Eq. (22d). Actually, we have found that the longitudinal paiTof the PN 
displacement vector first appears of the same order of its magnitude. clP Thus the 
existence of the growing transverse part in the PN displacement vector is interesting 
because such a mode is present in neither Newtonian theory nor the gauge invariant 
linear theory. 

Buchert has shown that the Newtonian trajectory field has no growing trans¬ 
verse mode up to third order in A when the initial velocity field is assumed to be 
irrotational.0 In his work, the solution of the transverse part at order 

A^ takes the form 


9 


T 

22 : 2 : 


p\3ri/t\2 3/t\4/3 ip\2/3 ^ 4 /^ ^-1/3 4 ^^-!- 

(40a) 


with 




(iV):7(3) _ ^(1) ^(2) 


F 


-F 


(1) ,^,(2) 


(1) .r,(2) 


+ F^ F 


-F 


( 1 ) 1 ^ 7 ( 2 ) 


Nlijk'^Nljk "^NlkjjNlki ^ ^N\ki N\kjj N\kj N\kij 


(40b) 


He has pointed out that the existence of the above transverse part ensures the con¬ 
servation of the vorticity flow along the fluid flow up to order A^. 

To see the importance of the PN effect on Newtonian dynamics, let us make 

o 

a simple order estimation. Consider the initial density fluctuation Sg{i) with the 
characteristic length I in the comoving frame. Note that due to o(to), we may regard 
the length I as the physical scale at the present time to- By using Eqs. (22d), (22e), 
(38), (39),and (^), the above transverse parts of the Newtonian displacement vector 
and the PN displacement vector may be estimated respectively roughly at some time 
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( 41 ) 

where dg(^i){tQ) is defined by (l + z/) 5g(i) representing the present value of the density 
contrast for the scale I extrapolated by the linear theory. This expression (41) 
suggests the following interpretations. As long as we consider the evolution of the 
transverse part with scale I much smaller than the present horizon scale cto, it is 
described well only by Newtonian theory, including the transverse part, because the 
Newtonian transverse part remains much larger than the PN transverse 

part from the initial time. This may not be correct for fluctuations 

with larger scale. In fact, from Eq. (41) we can evaluate the redshift zm when the 
magnitude of the Newtonian transverse part begins to become larger than that of 
the PN transverse part for the density fluctuation field with scale 1: 

1 + ZAT ~ (<^g{o(^o)) ^ ■ (42) 

Noting that cto is on the order of the present horizon scale and 2000/1 ^Mpc, it is 
found that this redshift is smaller for fluctuations with larger scale. This means that 
the PN transverse part might play an important role for the early evolution of large- 
scale structure. Since some fluctuations with large scale become localized smaller in 
comparison with the horizon scale, which grows larger, it is naturally described more 
accurately by Newtonian dynamics later. For example, SDSS is expected to survey a 
cubic region of several hundred megaparsecs length. Then the above estimate gives 
zat ~ 4 for fluctuations with scale 50/i“^Mpc at the present time, if we estimate 
the magnitude of the present density contrast as (5^(50j(/o) ~ 10“^, according to the 
power spectrum of the density contrast in the standard CDM modeli^ normalized 
by the COBE observation. 

We can solve for the trajectory field, namely the integral curve of velocity field 
V in our coordinates, up to the 1st PN order, so it will be interesting to see how 
the relativistic effect on Newtonian dynamics appears as a characteristic pattern in 
the large-scale structure formation. It will be easy to investigate this problem ex¬ 
panding the Newtonian simulation based on the Lagrangian perturbation approach. 
Also, as stated above, the PN characteristic effect might appear on larger scales. In 
particular, the effecimnay induce secondary temperature anisotropies of CMB, as the 
Rees-Sciama effect, E2P and potentially be observable. This reason in the anisotropy 
may be interpreted as follows. As stated in §2, we foliate the spacetime by simul¬ 
taneous surfaces of the basic observers who see no dipole component of CMB using 
the (3-|-l) formalism. Thus, it would be natural for the existence of the shift vector 
to influence the anisotropy of CMB. These works are now in progress. 

Since the observable region of the large-scale structure increases steadily, we 
expect that our work will play some important role in the future. 
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